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Background: Estrogen is known to decrease the risk of colon cancer in postmenopausal women, 
and may exert its actions by decreasing interleukin-6 (IL6) production via stabilization of the 
transcription factor nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB). 
Estrogens are biosynthesized by CYP19A1 (aromatase), so it is possible that genetic variations 
in CYP19A1 influences the risk of colon cancer by altering expression of CYP19A1. Further, 
studies on gene-gene interactions suggest that single nucleotide polymorphisms in one gene 
may affect expression of other genes. The current study aims to explore the role of CYP19A1 
single nucleotide polymorphisms on CYP19A1, NFkBI and IL6 gene expression. 
Methods: Phenotype-genotype associations, cross-associations between genes, and haplotype 
analyses were performed in both normal human colon (n=82) and liver (n=238) samples. 
Results: CYP19A1 rsl0459592, rsl961177, and rs6493497 were associated with CYP19A1 
expression in colon samples (P=0.042, P=0.041, andP=0.013, respectively). CYP19A1 single 
nucleotide polymorphisms (rsl2908960, rs730154, rs8025191, and rsl7523880) were cor- 
related with NFkBI expression (P=0.047, P=0.04, P=0.05, and P=0.03, respectively), and 
CYP19A1 rsl 1856927, rs2470152, and rs2470144 (P=0.049, P=0.025, P=0.047, respectively) 
were associated with IL6 expression in the colon. While rs730154 and rsl7523880 could not 
be analyzed in the liver samples, none of the other associations with the colon were replicated 
in the liver samples. Haplotype analysis revealed three separate haplotypes of the CYP19A1 
single nucleotide polymorphism that were significantly associated with CYP19A1, NFkBI, 
and IL6 gene expression. 

Conclusion: CYP19A1 single nucleotide polymorphisms are associated not only with CYP19A1 
expression but also with NFkBI and IL6 expression. These data demonstrate the possible 
functional consequences of genetic variation within the CYP19A1 gene on other genes in a 
biologically plausible pathway. 
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Introduction 

Colorectal cancer is the third most common cancer and cause of cancer death in women.' 
The link between estrogen and the risk of colon cancer has been well established.^"* 
Hormone replacement therapy has been shown to decrease the risk of colon cancer 
in postmenopausal women,' with two separate meta-analyses demonstrating a 20% 
decrease in risk of colon cancer in ever users. ''•^ 

One mechanism by which estrogen exerts its protective action may be via modulation 
of the immune system response.** A recent study in mice indicated that estrogen suppresses 
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inflammation-associated colon cancer, even after initiation of 
colonic DNA damage by administration of azoxymethane 
and dextran sulfate sodium.' Specifically, estradiol negatively 
regulates the inflammatory cytokine, interleukin-6 (IL6)."' " 
It has been demonstrated that IL6 levels increase as Dukes' 
stage increases inpatients with colorectal cancer, so increased 
IL6 expression may play a role in the risk of colon cancer.'"^ 
CYP19A1 (aromatase), which facilitates the biosynthesis of 
estrogens, is expressed by colorectal adenocarcinoma cell 
lines," and colon carcinoma has increased expression of 
CYP19A1 mRNA relative to non-neoplastic colon mucosa.''' 
While the mechanisms of this difference in expression have 
not been identified, it has been postulated that an inflamma- 
tory response process may be involved, with expression of 
IL6 and C YP 1 9 A 1 activity having been positively correlated 
in macrophage-rich tissue" and IL6 and other inflammatory 
cytokines having been shown to increase the expression of 
CYP19A1. This increases local estrogen levels,"" which can 
inhibit binding of nuclear factor kappa-light-chain-enhancer 
of activated B cells (NFkB) to inflammatory gene promoters. " 
NFkB is a transcription factor that has binding sites in or near 
many cytokine promoter regions,'^ and is closely associated 
with the development of cancer, angiogenesis, invasion, and 
metastasis." NFkB is known to mediate the synthesis of IL6,^'' 
indicating a possible feedback loop. 

Single nucleotide polymorphisms (SNPs) in CYPI9A1 
have been linked to a risk of colon cancer.'' Interactions 
between SNPs in pairs of genes, such as CYP19A1 and 
NFkBI, have also been shown to affect the risk of colon can- 
cer.^^ The functionality of these associations, as determined by 
gene expression, has not been explored. The primary purpose 
of this study was to validate findings from a large case-control 
study^^ by exploring the gene expression functionality of (and 
possible interactions between) SNPs in CYP19A1 and gene 
expression of CYP19A1, NFkBI, and IL6 in normal colon 
samples. Haplotype analysis of multiple CYP19A1 SNPs was 
performed to understand better the broader variation in the 
CYP19A1 gene and how it may influence expression. 

Materials and methods 

Tissue samples 

Deidentified normal frozen colon and liver tissue samples 
were obtained from the Cooperative Human Tissue Network 
and stored at -80°C. For colon tissues (n=82), 54% were from 
male patients and 46% were from female patients. The mean 
patient age was 60.48 (range 17-92) years, and the samples 
were from patients of Caucasian (n=5 1), African American 
(n=23), Asian (n=l), or unknown (n=7) origin. 



For the liver samples (n=238), the mean patient age was 
58.79 (range 1-102) years, and the samples were from patients 
of Caucasian (n=195), African American (n=21), Hispanic 
(n=2), or other/unknown (n=20) origin. Fifty-three percent of 
the liver samples came from male patients, 45%o from female 
patients, and 2% from patients of unknown sex. 

Reverse transcription 

and quantitative real-time PGR 

To increase the yield in normal colon tissue samples, total 
RNA was isolated utilizing Trizol reagent (Invitrogen, Grand 
Island, NY, USA) for homogenization, and the RNEasy Mini 
kit (Qiagen, Valencia, CA, USA) for isolation using a proto- 
col developed by Mauricio Rodriquez-Lanetty (unpublished) 
but with minor alterations. Briefly, tissue samples (~25 mg) 
were homogenized in 1 50 |J,L of Trizol using a Bullet Blender 
(Next Advance, Averill Park, NY, USA) and stainless steel 
beads. The homogenate was placed in a new vial with 450 |J,L 
of Trizol. After adding 100 jlL of chloroform, the vials were 
shaken, incubated for 2 minutes at room temperature, and 
centrifuged, after which the supernatant was placed in a new 
vial. An equal portion of 100% ethanol was added, and the 
mixture was placed in an RNEasy spin column. The RNA was 
washed and eluted according to the RNEasy protocol. Total 
DNA was isolated from normal colon samples and total RNA 
and DNA were isolated from normal liver tissue samples 
using the AUPrep DNA/RNA/Protein mini kit (Qiagen). 

First strand cDNA synthesis was performed using a High 
Capacity RNA-to-cDNA Kit (ABI, Carlsbad, CA, USA) on 
500 ng of total RNA as measured by an RNA 6000 Nano kit 
(Agilent, Santa Clara, CA, USA). Quantitative real-time PCR 
reactions were performed on a 7900HT Fast Real-Time PCR 
System (ABI) using Taqman primer/probe sets and Taq- 
man Fast Universal PCR Master Mix no AmpErase* UNG 
(ABI). Experiments were run as per the manufacturer's 
protocol (except as noted) in quadruplicate on CYP19A1 
colon cDNA and in triplicate on NFkBI and IL6 cDNA 
and CYP19A1 liver cDNA diluted 1:10 for 50 PCR cycles. 
Samples were normalized to p-actin, discarding those with 
P-actin cycle thresholds (Cts) that were >30 (n=12 for colon 
CYP 1 9A 1 analysis, n=5 for colon IL6 analysis, n=4 for colon 
NFkBI analysis; and n=8 for all liver analyses). Gene of 
interest Cts that were &40 or that were undetermined were set 
to 40 (n=13 for colon CYP19A1 analysis, n=2 for colon IL6 
analysis, n=86 for liver CYP19A1 analysis, and n=9 for liver 
IL6 analysis). Two observations were removed from the liver 
CYP19A1 analysis because the standard deviation among the 
CYP19A1 triplicates was >1; one observation was similarly 
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removed from the liver IL6 analysis. P-actin was used as the 
housekeeping gene because it has been shown that structural 
housekeeping genes such as p-actin have less variation in 
normal colon and liver tissues than metabolic housekeeping 
genes such as glyceraldehyde 3-phosphate dehydrogenase.^'' 

SNP genotyping of colon samples 

The process of choosing TagSNPs for analysis has been 
described previously.-^ All markers were genotyped using a 
multiplexed bead array assay format based on GoldenGate 
chemistry (lUumina, San Diego, CA, USA). A genotyping call 
rate of 99.93% was achieved. Blinded internal replicates repre- 
sented 1 .6% of the sample set. The duplicate concordance rate 
was 99.996%. In this study, 25 SNPs from CYP19A1, eleven 
from NFkBI, and five from IL6 were examined. 

SNP genotyping of liver samples 

Genotyping for SNPs found to be significant in the colon 
was performed using Taqman SNP Genotyping Assays (Life 
Technologies, Carlsbad, CA, USA) on the 7900HT Fast Real- 
Time PCR System (ABI). The PCR protocol was adapted as 
follows: AmpliTaq activation (95°) for 10 minutes, followed 
by 50 cycles of denature (95°) for 1 5 seconds and anneal (58°) 
for 90 seconds. Samples were automatically called with the 
post-PCR allelic discrimination protocol, with a genotyping 
success rate of >95%. One SNP (rs730154) did not have 
a Taqman primer available, so was not analyzed in liver 
samples. CYP19A1 rsl7523880 was not in Hardy-Weinberg 
equilibrium, so was not analyzed either. 

Statistical analysis 

Effects of individual SNPs on expression 

Statistical analyses were performed using SAS version 
9.3 software (SAS Institute, Cary, NC, USA). Tests for 
Hardy-Weinberg equilibrium and linkage disequilibrium 
measures were calculated and stratified by race using the 
ALLELE procedure. Median 2'^ACt (ACt = control Ct - gene 
of interest Ct)^^ values were calculated by genotype initially 
assuming a codominant model, and the best fitting inheritance 
model is presented. P- values based on nonparametric analysis 
of variance tests (Wilcoxon rank-sum and Kruskal-Wallis 
rank-sum tests) were used to detect differences in distribution 
of expression levels by genotype. 

Haplotypes 

CYP19A1 SNPs that appeared to be significantly associ- 
ated (at a=0.05) with expression were used to construct 
haplotypes. Haplotypes were not created based on linkage 



disequilibrium bins as tagSNP markers were chosen from 
different linkage disequilibrium bins. Haplotypes were 
created based on results from colon samples as follows: 
all CYP19A1 SNPs that were significantly associated with 
CYPI9AI expression (n=3); all CYPI9A1 SNPs that were 
significantly associated with NFkBI expression (n=4); and 
all CYP19A1 SNPs that were significantly associated with 
IL6 expression (n=3). 

The haplotype reconstruction software, PHASE,^*" was 
used to predict the population frequencies and the "best pair" 
of haplotypes for each individual. This program provides 
a "best pair" of haplotypes that most likely explains the 
genotypes of that subject with the corresponding probability 
of having that particular haplotype for each individual. If 
haplotypes for each individual are uncertain, then the pro- 
gram assigns alternative pairs of haplotypes with respective 
probabilities. PHASE output was read into SAS version 9.3 
which created one variable for each haplotype with three 
possible values, ie, 0 (no copies of that haplotype), 1 (one 
copy of that haplotype), or 2 (two copies of that haplotype). 
Association of each haplotype with gene expression was 
examined again using nonparametric analysis of variance. 
A P-value of less than 0.05 (two-sided) was considered to 
be statistically significant. Any haplotypes found to be sta- 
tistically significant in normal colon samples were tested in 
normal liver samples as well. 

Results 

Effects of CyP/9A/ SNPs 
on CYP 1 9A I expression 

Three CYP19A1 SNPs were significantly associated with 
changes in gene expression in colon samples (Table 1). The 



Table I Association of CYP 1 9A I expression and CYP/9A/ SNPs 
in colon samples 



SNPs 


Frequency 
(n) 


Median 
(l-^ACt) 


Fold 
change 


P-value* 


rs 10459592 










TT/TG 


53 


9.95E-05 


4.7 


0.042 


GG 


17 


4.65 E-04 






rsl96l 177 










CC 


5! 


I.53E-04 


1.8 


0.040 


CT/TT 


19 


8.37E-05 






rs6493497 










GG 


53 


I.53E-04 


1.9 


0.013 


GA/AA 


17 


8.I0E-05 







Notes: rsl96l 177 and rs6493497 r^l among Caucasians and 0.58 among African- 
Americans. ^Nonparametric analysis of variance {Kruskal-Wallis or Wilcoxon rank 
sum test). 

Abbreviation: SNPs, single nucleotide polymorphisms. 
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CYP19A1 rsl0459592 homozygous variant genotype (GG) 
was associated with 4.7-fold higher expression (/'=0.042) of 
CYP19A1 compared with the heterozygous or homozygous 
common (TT/TG) genotype. The rs 1 96 1 1 77 CT/TT genotype 
was associated with 1 .8-fold lower gene expression (P=0.040) 
as compared with the CC genotype. Having one or more cop- 
ies of the rs6493497 variant allele (GA/AA) was associated 
with 1.9-foId lower (/'=0.013) CYP19A1 expression. No 
significant associations between CYP19A1 SNPs and mRNA 
expression levels were found in the liver samples. 

Cross-associations between 
aPl9Al SNPs and NF/cfi/ and /L6 
expression in colon samples 

The CYP19A1 rsl2908960 homozygous variant genotype 
(AA) was associated with 2.0-fold lower (/'=0.047) NFkBI 
expression as compared with the heterozygous or homozy- 
gous common (GG/GA) genotype. Having one or more 
copies of the CYP19A1 rs730154 or rs8025191 variant allele 
(AG/GG) was associated with 1 .6-fold higher (_P=0.036) and 
1.7-fold lower (although this only approached significance 
at P=Q.Q5\) NFkBI expression, respectively, as compared 
with the homozygous common (AA) genotype. Having one 
or more copies of the CYP19A1 rsl7523880 variant allele 
(CA/AA) was associated with 1.8-fold lower (P=0.031) 
NFkB 1 expression as compared with the homozygous com- 
mon (CC) genotype (Table 2, A). 

The CYP19A1 rsl 1856927 homozygous variant genotype 
(GG) was associated with 3.1-fold higher (P=0.049) IL6 
expression compared with the heterozygous or homozygous 
common genotype (TT/TG). Having two of the variant alleles 
in CYP19A1 rs2470152 (TT) was associated with 3.7-fold 
higher (/'=0.025) IL6 expression compared with the heterozy- 
gous or homozygous common genotype (CC/CT), while hav- 
ing one or more of the variant allele in CYP19AI rs2470144 
(GA/AA) was associated with 3.4-fold higher {P=0M7) 
IL6 expression (Table 2, B). A summary of /"-values for all 
single nucleotide polymorphisms explored in colon samples 
can be found in Table SI. 

Cross-associations between CYPI9AI 
SNPs and IL6 expression in liver samples 

One or two copies of the variant (A) allele of CYP19A1 
rsl 2908960 or rs2470144 was significantly associated with 
3.5-fold and 2.2-fold lower IL6 expression compared with 
the homozygous common (GG) genotype (^=0.004 and 
P=0.022, respectively. Table 2, C). A summary of P-values 



Table 1 Cross-associations of CYPI9AI genotypes with IL6 and 


NFkBI gene expression in 


colon and liver 


samples 




CYPt9A\ 


Frequency 


Median 


rolo 


P-value* 


SNPs 


(n) 


(l-^ACt) 


change 




(A) NFkBI expression in colon 






rs 1 2908960 










GG/GA 


65 


0.079 


2.0 


0.047 


AA 


13 


0.040 






rs730IS4 










AA 


39 


0.052 


1.6 


0.036 


AG/GG 


39 


0.085 






rs8025 1 9 1 










AA 


5! 


0.085 


1.7 


0.051 


AG/GG 


27 


0.050 






rs 1 7523880 










CC 


5! 


0.082 


1.8 


0.031 


CA/AA 


27 


0.045 






(B) IL6 exp 


resslon in colon 






rsl 1856927 










TT/TG 


6! 


0.003 1 


3.1 


0.049 


GG 


16 


0.0094 






rs2470l52 










CC/CT 


65 


0.0033 


3.7 


0.025 


TT 


12 


0.0121 






rs2470l44 










GG 


33 


0.0023 


3.4 


0.047 


GA/AA 


44 


0.0079 






(C) IL6 expression In liver 








rsl 2908960 










GG 


67 


0.0149 


3.5 


0.004 


GA/AA 


106 


0.0043 






rs2470l44 










GG 


66 


0.0145 


2.2 


0.022 


GA/AA 


139 


0.0066 







Note: *Nonparametric analysis of variance (Kruskal-Wallis or Wilcoxon rank sums 
test). 

Abbreviations: IL6, interteukin-6; SNPs, single nucleotide polymorphisms: NFkB, 
nuclear factor kappa-light-chain-enhancer of activated B-cells. 



for all single nucleotide polymorphisms explored in liver 
samples can be found in Table S2. 

Haplotype analysis 

The CYP19A1 rsl0459592 T>G, rsl961177 C>T, and 
rs6493497 G>A haplotype TTA (one copy) was sig- 
nificantly associated with CYP19A1 expression (P=0.003, 
frequency 0.20; Table 3, A). Haplotype construction of 
CYP19A1 rsl2908960 G>A, rs730154 A>G, rsl7523880 
OA, and rs8025191 A>G demonstrated that two copies 
of GGCA (frequency 0.12) was associated with increased 
NFkBI expression (/'=0.049; Table 3, B). A single copy of 
another CYP19A1 haplotype (rsl 1856927 T>G, rs2470152 
OT, rs2470144 G>A) was associated with IL6 expres- 
sion (TCG, P=0.011, frequency 0.29), with a different 
combination of GCG approaching significance at P=0.053 
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Table 3 Association of haplotypes of CYPI9A I SNPs with selected 
gene expression in colon 



CYPI9AI 
SNPs 



Frequency 



Median 
(2'^ACt) 



SD 



P-value* 



(A) CYP 1 9A I expression 

rs 1 0459592 T>G 
rsl96l 177 OT 
rs6493497 G>A 
Haplotype 
TTA 

0 copyt 0.80 

1 copy 0.20 

(B) NFkBI expression 

rs 1 2908960 OA 
rs730l54 A>G 
rs 1 7523880 OA 
rs8025l9l A>G 
Haplotype 
GGCA 



1 .8E-4 
4.8E-5 



0.009 
0.0002 



0.003 



0 copy 


0.51 


0.048 


0.099 


0.049 


1 copy 


0.37 


0.085 


0.134 




2 copies 


0.12 


0.114 


0.126 





(C) IL6 expression 

rsl 1856927 T>G 
rs2470l52 OT 
rs2470l44 G>A 
Haplotypes 
TCG 



0 copy 


0.71 


0.0078 


0.151 


1 copy 


0.29 


0.0019 


0.018 


GCG 








0 copy 


0.78 


0.0047 


0.123 


1 copy 


0.16 


0.0014 


0.010 


2 copies 


0.06 


0.067 


0.273 



0.01 I 



0.053 



Notes: Haplotypes created based on results from colon samples; *Nonparametric 
analysts of variance (Kruskai-Waliis or Wilcoxon rank sum test); ^0 copy, no copies 
of the haplotype: I copy, one copy of the haplotype; 2 copies, two copies of the 
haplotype. 

Abbreviations: IL6, interleukln-6; SNPs, single nucleotide polymorphisms; NFkB, 
nuclear factor kappa-light-chain-enhancer of activated B-cells; SD, standard deviation. 



(Table 3, C). A summary of haplotype analyses can be found 
in Table S3. 

Discussion 

Understanding the interplay between estrogen synthesis 
genes and inflammatory genes may help to explain the 
mechanisms behind the role of estrogen in reducing the 
risk of colon cancer. The estrogen synthesizer CYP19A1 
(aromatase, coded by the CYP19A1 gene), the transcrip- 
tion factor NFkB (coded by NFkBI), and the inflammatory 
cytokine IL6 (coded by IL6) may work together to create a 
feedback loop for homeostasis. An increase in CYP19A1 
would lead to an increase in estrogen, which could lead to 
decreased IL6 production by stabilization of NFkB. The 



purpose of this study was to explore the gene expression 
functionality of, and possible interactions between, SNPs in 
CYP19A1 and gene expression of CYP19A1, NFkBI, and 
IL6 in normal colon samples. 

Variants in CYP19A1 (rsl0459592, rsl961177, 
rs6493497) were significantly associated with CYP19A1 
gene expression in colon samples. CYP19A1 rsl961177 
has previously been associated with the risk of colon cancer 
and rectal cancer. Transcription factors and other regula- 
tory molecules vary in expression across organs and tissues, 
resulting in tissue-specific regulation. Since the liver is the 
site of first-pass metabolism for many colon carcinogens, 
genotype-phenotype correlations in DNA and mRNA from 
238 normal liver samples were also examined. In contrast 
with colon tissue, none of the CYP19A1 SNPs that had a sig- 
nificant association in the colon were significantly associated 
with CYP19A1 expression in liver samples. The differences 
seen for CYP19A1 SNP associations with CYP19A1 gene 
expression between the organs could be due to differences in 
expression of transcription factors and regulatory molecules 
between colon and liver tissues. 

When cross-associations were examined in the colon, 
CYP19A 1 rs 12908960, rs730 1 54, rs8025 1 9 1 , and rsl 7523880 
were associated with NFkBI expression. In addition, 
CYP19A1 rsl 1856927, rs2470152, and rs2470144 were 
significantly associated with IL6 gene expression. While no 
single SNP in CYP19A1 influenced expression in all three 
genes, these data indicate there may be some cross-talk 
between genes that affect function. It is possible that the 
cross-gene associations identified could occur by chance, 
and rigorous mechanistic testing needs to be performed to 
verify these observations. 

CYP19A1 SNPs found to be significant in the colon were 
also explored for cross-association in the liver. CYP19A1 
rsl2908960 and rs2470144 were associated with IL6 
expression in the liver samples. While one or more copies 
of the variant allele genotype (GA/AA) for rs2470144 were 
associated with increased expression in colon samples, the 
variant allele was associated with decreased expression in 
liver samples. Given this difference in association and the fact 
that no other cross-associations were found to be significant 
in the liver, these associations could be organ-specific. 

Each of the genes is located on a separate chromosome 
(CYP19A1 on 15, NFkBI on 4, and IL6 on 7). Haplotypes 
were developed with CYP19A1 SNPs, and expression levels 
of CYP19A1, NFkBI, and IL6 mRNA were explored in 
colon samples. Haplotype/activity analysis demonstrated 
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that haplotype TTA constructed from CYP19A1 rsl0459592 
T>G, rs 1961 177 C>T, andrs6493497 G>A was associated 
with CYP19A1 mRNA expression. The CYP19A1 haplotype 
GGCA (two copies) constructed from rsl2908960 G>A, 
rs730154 A>G, rsl7523880 OA, and rs8025191 A>G 
demonstrated a significant influence on expression of NFkB 1 . 
The haplotype TCG constructed from CYP19A1 rs 1 1 856927 
T>G, rs2470152 C>T, andrs2470144 G>A was associated 
with IL6 mRNA expression. The haplotype GCG (con- 
structed from the same SNPs) also showed an association 
with IL6 mRNA expression. This association, however, only 
approached statistical significance (/'=0.053). These data 
suggest that CYP19A1 SNPs, as well as haplotypes, affect 
NFkBI and IL6 expression, as well as CYP19A1 expression. 
The CYP19A1 haplotype rsl2591359, rsl7523880, 
rsl961 177, and rs3751591 was reported by Slattery et al to 
be associated with risk of colon cancer (global P=0.0059). 
This haplotype was tested in normal colon samples in our 
study, and one construct, ACTT, was significant (/'=0.008, 
frequency 0.06). This finding supports the association 
reported by Slattery et al.-^ 

Further studies of the coregulation of these genes and 
their effect on risk of colorectal cancer need to be performed 
in a large population sample to confirm that these are not 
false findings based on the small number of colon samples. 
Samples should also include individuals from different ethnic 
backgrounds. Other mechanisms also need to be considered, 
such as the role of transcription factors in these associations. 
Further study of how NFkBI may affect the associations 
shown here is warranted, as well as exploration of other 
suitable transcription factors. For example, the transcription 
factor, NF-IL6, has been associated with both CYP19AP' 
and (synergistically withNFKB) IL6^'* expression in endome- 
trial stromal cells and murine PI 9 cells, respectively. These 
associations have not, however, been tested in colon samples. 
Future studies will include testing this transcription factor at 
the genetic level, as well as examining the combined effect 
of these genes and their SNPs on the risk of colon cancer in 
case-control studies. 

These data demonstrate the functional consequences of 
genetic variation within one gene on interconnected genes 
within a biologically plausible pathway. They also indicate 
that there are cross-associations between CYP19A1 SNPs and 
NFkB 1 and IL6 expression in colon samples, although these 
findings were not validated in liver samples. Further research 
may help to elucidate the mechanisms behind the relationships 
described here. Understanding the estrogen-metabolizing 
CYPl 9 A 1 gene and its interactions with other pertinent genes 



could lead to a better understanding of colon cancer and to 
more effective prevention and treatment options. 
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Table SI Association of CYPI9A1, NFkBI, and IL6 expression 
and CYPI9AI, SNPs in colon samples 

SNPs Location CYPI9AI NFkBI IL6 

expression expression expression 
P-value P-value P-value 
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Abbreviations; IL6. interleukin-6; SNPs, single nucleotide polymorphisms; NFkB, 
nuclear factor kappa-light-chain-enhancer of activated B-cells. 



Table S2 Association of CYPI9AI, NFkBI, and IL6 expression 
and CYPI9AI SNPs in liver samples 

SNPs CYPI9AI NFkBI iL6 

expression expression expression 
P-value P-value P-value 
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Abbreviations: SNPs, single nucleotide polymorphisms; IL6, lnterleukin-6; NFkB, 
nuclear factor kappa-light-chain-enhancer of activated B-cells. 
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Table S3 Association of haplotypes of CYPI9AI SNPs with Table S3 (Cont/nued) 
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